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Abstract 

We consider an invariant quantum Hamiltonian H = —Alb + V in the L 
space based on a Riemannian manifold M with a discrete symmetry group T. 
Typically, M is the universal covering space of a multiply connected manifold 
M and T is the fundamental group of M. To any unitary representation A of 
r one can relate another operator on M = M/T, called H\, which formally 
corresponds to the same differential operator as H but which is determined 
by quasi-periodic boundary conditions. We give a brief review of the Bloch 
I/"") I decomposition of H and of a formula relating the propagators associated to the 

Hamiltonians H\ and H . Then we concentrate on the example of the Aharonov- 
Bohm effect with two vortices. We explain in detail the construction of the 
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propagator in this case and indicate all essential intermediate steps. 
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! 1 Introduction 

w 

Suppose that there is given a connected Riemannian manifold M with a discrete sym- 
metry group T. Let us consider a T-periodic Hamilton operator in L 2 (M) of the form 
H = — Alb + V where Alb is the Laplace-Beltrami operator and V is a T- invariant 
bounded real function on M. To any unitary representation A of T one can relate 
another operator on M = M/T, called H\, which formally corresponds to the same 
differential operator as H but which is determined by quasi-periodic boundary condi- 
tions. In the framework of the Feynman path integral there was derived a remarkable 
formula relating the propagators /Cf (x, x ) and JC t (x, x ) associated respectively to the 
Hamiltonians if a and H [TJ |2] . An analogous formula is also known for heat kernels 
There exists also an opposite point of view when one decomposes the operator H 
into a direct integral with components H\ where A runs over all irreducible unitary 
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representations of T [H EJ E] . The evolution operator then decomposes correspond- 
ingly. This type of decomposition is an essential step in the Bloch analysis. Let us also 
note that an alternative approach to the Bloch analysis, based on a more algebraic 
point of view, has been proposed recently in [TJ. 

The both relations, the propagator formula on the one hand and the generalized 
Bloch decomposition on the other hand, are in a sense mutually inverse [S]. In the 
current paper we give a brief review of basic results concerning this relationship. The 
main purpose of this contribution is, however, a detailed discussion of an application 
of the formula for propagators. We consider the example of the Aharonov-Bohm effect 
with two vortices. In this case M is identified with the universal covering space of the 
plane with two excluded points and T is the fundamental group of the same manifold. 
This problem has been already treated by one of the authors quite a long time ago in 
[9]. But even in more recent papers one can encounter a discussion of the problem itself 
[TO] as well as of the involved methods [UJ, [12]. It should be stressed that article [9] 
is a very brief letter in which only the final formula is presented without any detailed 
hints about its derivation. This is so even though some steps in the derivation are 
rather complicated and in no way obvious. Since these details were never published 
anywhere and since without them the resulting formula may look a bit obscure we aim 
here to fill in this gap by explaining the approach more carefully and by indicating 
the necessary intermediate steps. 

The paper is organized as follows. In Section [5] we give a brief review of basic 
results concerning the relationship between the generalized Bloch analysis and the 
formula for propagators associated to periodic Hamiltonians. In Section |3] we explain 
the construction of the propagator on the universal covering space in the case of the 
Aharonov-Bohm effect with two vortices. In Section H] we discuss the application of 
the propagator formula in this particular case. 

2 Propagators associated to periodic Hamiltonians 
2.1 Periodic Hamiltonians 

Let M be a connected Riemannian manifold with a discrete and at most countable 
symmetry group T. The action of V on M is assumed to be smooth, free and proper 
(also called properly discontinuous). Denote by jl the measure on M induced by the 
Riemannian metric. The quotient M = M /V is a connected Riemannian manifold 
with an induced measure fi. This way one gets a principal fiber bundle tt : M — *■ M 
with the structure group T. The L 2 spaces on the manifolds M and M are everywhere 
tacitly understood with the measures /i and jl, respectively. 

Typically, M is the universal covering space of M and F = 7Ti(M) is the fun- 
damental group of M. For example, this is the case when one is considering the 
Aharonov-Bohm effect. 

To a unitary representation A of T in a separable Hilbert space one relates the 
Hilbert space formed by A-equivariant vector- valued functions on M. This means 
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that any function ip G M 3 ^ is measurable with values in and satisfies 

Vs G T, ip(s ■ y) = A(s)ijj(y) almost everywhere on M. 

Moreover, the norm of if) induced by the scalar product is finite. If ipi,ip2 £ then 
the function y \— > (ipi(y), faiy)) defined on M is T-invariant and so it projects to a 
function defined on M, and the scalar product is defined by 

(^1,^2)= / s^ 3 (x)dfi(x). 

J M 

As already announced, our discussion concerns T-periodic Hamiltonians on M of 
the form H = —Alb + V where Alb is the Laplace-Beltrami operator and V(y) is a 
T-invariant measurable bounded real function on M. Here we accept the Friedrichs 
extension as the preferred self-adjoint extension of semibounded symmetric operators 
defined on test functions. 

To the same differential operator, — Alb + V, one can relate a selfadjoint operator 
if a in the space J^f\ for any unitary representation A of T. Let us define <E>a : Cq°(M)® 
Sf A -> JZk by 

G C °°(M), G J^ A , ($ A ^ <g> v) (y) = • V) A(s~>- 

ser 

Since the action of T is proper, the vector- valued function $a f ® f is smooth. More- 
over, $a V 9 ® u is A-equivariant and the norm of $a ® v in ^a. is finite. Furthermore, 
the range of $a is dense in J#a. The Laplace-Beltrami operator is well defined on 
Ran($A) and it holds 

Alb^kW ® «] = $a[A lb ^ ® «]. 

One can also verify that the differential operator —Alb is positive on the domain 
Ran($A) C J^a- Since the function V(y) is T-invariant, the multiplication operator 
by V is well defined in the Hilbert space The Hamiltonian H\ is defined as the 
Friedrichs extension of the differential operator —Alb + V considered on the domain 
Ran$ A . 

2.2 A generalization of the Bloch analysis 

Let f be the dual space to V (the quotient space of the space of irreducible unitary 
representations of V). In the first step of the generalized Bloch analysis one decomposes 
H into a direct integral over T with the components being equal to if a- As a corollary 
one obtains a similar relationship for the evolution operators U(t) = exp(-itH) and 
UA.(t) = exp(-itHx), tel. To achieve this goal a well defined harmonic analysis on 
the group T is necessary. 

It is known that the harmonic analysis is well established for locally compact groups 
of type I [TJ]. So all formulas presented bellow are perfectly well defined provided T 
is a type I group. A countable discrete group is type I, however, if and only if it has 
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an Abelian normal subgroup of finite index [HJ Satz 6]. This means that there exist 
multiply connected configuration spaces of interest whose fundamental groups are not 
of type I. For example, the fundamental group in the case of the Aharonov-Bohm 
effect with two vortices is the free group with two generators and it is not of type I. 
Fortunately, in this case, too, there exists a well defined harmonic analysis |15j . 

Let us recall the basic properties of the harmonic analysis on discrete type I groups 
|13j . In that case the Haar measure on T is chosen as the counting measure. Let dm 
be the Plancherel measure on T. Denote by ^(-^a) = ® the Hilbert space 
formed by Hilbert-Schmidt operators on JzfA (=^a is the dual space to Jz?a)- The 
Fourier transformation is defined as a unitary mapping 

For / G L 1 ^) C L 2 (r) one has 

^[/](A) = ^/( S )A( S ). 
ser 

Conversely, if / is of the form / = g * h (the convolution) where g,h G ^(T), and 
/ = &[f\ then 

/( S ) = ^Tr[A(s)*/(A)]dm(A). 

It is known that if T is a countable discrete group of type I then dimJzfA is a 
bounded function of A on the dual space T [HJ Korollar I]. Using the unitarity of the 
Fourier transform one finds that 

m(t) < J dim«jSf A dm(A) = 1. 

The following rule satisfied by the Fourier transformation is also of crucial importance: 
V. G T, V/ G L 2 (r), &[f(s ■ g)](A) = A( S - 1 )^[/(( 7 )](A). 
Now we are going to construct a unitary mapping 

r® 

$ : L 2 (M) -> J J^ A * ® ^ A dm(A) 

which makes it possible to decompose the Hamiltonian H. Observe that the tensor 
product J^a^^^a can be naturally identified with the Hilbert space of l®A-equivariant 
operator- valued functions on M with values in Jzf^ ® Jz?a = -^(-^a)- For / G L 2 (M) 
and y G M set 

\/seT,f y (s) = f(s- 1 -y). 

The norm \\f y \\ in L 2 (T) is a T-invariant function of y G M, and the projection of this 
function onto M can be checked to be square integrable. Hence for almost all x G M 
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and all y G 7r _1 ({x}) it holds f y G L 2 (T). We define the component $[/](A), A G f, 
by the prescription 

In particular, if / G L l {M) n L 2 (M) then 

Equivalently one can define $ in the following way. For ip G C£°(M), v G Jz?a and 
2/ G M set 

$fo>](A)(y)i; = ($a¥>® (1) 

Then $ introduced in ([TJ) is an isometry and extends unambiguously to a unitary 
mapping. 

Finally one can verify the formula 

= J 1 <g> H A dm(A) 

which represents the sought Bloch decomposition. As a corollary we have 

re 

QUit)®- 1 = J l®£7 A (t)dm(A). (2) 

2.3 A construction for propagators associated to periodic 
Hamiltonians 

In equality ([2]), the evolution operator U(t) is expressed in terms of U\(t), A G T. 
It is possible to invert this relationship and to derive a formula for the propagator 
associated to H\ which is expressed in terms of the propagator associated to H. 

The propagators are regarded as distributions which are introduced as kernels of 
the corresponding evolution operators. Recall that by the Schwartz kernel theorem 
(see, for example, (TBI Theorem 5.2.1]), to every B G £$(L 2 (M)) there exists one and 
only one j3 G 3>'{M x M) such that 

Vpi,p 2 GC °°(M), P(<pi®<p 2 ) = (<p 1 ,B(p 2 ). 

Moreover, the map B \— > (3 is injective. One calls j3 the kernel of B. 

The kernel theorem can be extended to Hilbert spaces formed by A-equivariant 
vector-valued functions. In this case the kernels are operator-valued distributions. To 
every B G &(J#a) there exists one and only one (3 G @'{M x M) <g) ^(Jz?a) such that 

(Vi,P(<fi ® <f 2 )v 2 ) = ($ A ^1 <8> V!,B$ A (f 2 ® v 2 ). 

The distribution (3 is A-equivariant: 

Vs G T, j3{s ■ 2/ 2 ) = A(a)/%i, 2/2), s • 2/2) = 2/ 2 )A(s _1 ) 
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In this case, too, the map _B i — >• /? is injective. 

Denote by K t G &{M x M) the kernel of U(t) G ^(L 2 (M)), and by G 
3)'(M x M) ® £?(j£f A ) the kernel of £7 A (t) G 3§{3tf K ). Here and everywhere in this 
section, t is a real parameter. The kernel is A-equivariant: 

Ms G T, /C t A ( S ■ yi , |/ 2 ) = A(s)/C t V, 2/2), /Cf s • j/ 2 ) = /Cf ( m , y 2 )A( S - 1 ). 

First we rewrite the Bloch decomposition of the propagator (j2J) in terms of kernels. 
It is possible to prove that, for all <pi, <p 2 G Gq°(M), the function A h-> Tr[/Cf (^T <S> <£> 2 )] 
is integrable on T and 

£t(^i (8) ^2) = y Tr[/Cf ® ^2)] dm(A). 

An inverse relation was derived by Schulman in the framework of path integration 
[1, 2J and reads 

}Cf(x,y) = Y,Hs)}C t (s- 1 -x,y). (3) 

It is possible to give fl3J) the following rigorous interpretation. Suppose that y?i, y? 2 G 
Cq°(M) are fixed but otherwise arbitrary. Set 

F t (s) = /C t (^i(s -1 • y x ) ® ¥> 2 (ya)) for s G T, 

and 

G t (A) = /C t A (<^ <g> v? 2 ) G J^a) for A G f. 

One can show that F t G L 2 (T) and G t is bounded on f in the Hilbert-Schmidt norm. 
Recalling that m(f) < 1 we have \\G t (-)\\ G L x (f) n £ 2 (f ). In [8J it is verified that 

i^J^G*]. 

and, consequently, 

G t = &[F t \. (4) 
Rewriting (jlj) formally gives equality (J3j). 

3 The Aharonov-Bohm effect with two vortices: 
the propagator on the universal covering space 

3.1 A formula for the propagator 

The configuration space for the Aharonov-Bohm effect with two vortices is the plane 
with two excluded points, M = M 2 \{a, b}. This is a flat Riemannian manifold and the 
same is true for the universal covering space M. Let 7r : M — > M be the projection. It 
is convenient to complete the manifold M by a countable set of points A U B which lie 
on the border of M and project onto the excluded points, tt(A) = {a} and ir(B) = {b}. 
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M looks locally like M? but differs from the Euclidean space by some global features. 
First of all, not every two points from M can be connected by a geodesic segment. 
Fix a point x G M. The symbol D(x), as introduced below in (1151) . stands for the set 
of points y G M which can be connected with x by a geodesic segment. The domain 
D(x) is one sheet of the covering M — ► M. It can be identified with IR 2 cut along 
two halflines with the limit points a and b, respectively. Thus the border dD(x) is 
formed by four halflines. The universal covering space M can be imagined as a result 
of an infinite process of glueing together countably many copies of D(x) with each 
copy having four neighbors. 

The fundamental group of M, called T, is known to be the free group with two 
generators g a and g b . For the generator g a one can choose the homotopy class of a 
simple positively oriented loop winding once around the point a and leaving the point 
b in the exterior. Analogously one can choose the generator g^ by interchanging the 
role of a and h. One-dimensional unitary representations A of V are determined by 
two numbers a, /3, < a, (3 < 1, such that 

H9a) = e 2ma , A(g b ) = e 2 ^. 

The standard way to define the Aharonov-Bohm Hamiltonian with two vortices 
is to choose a vector potential A for which rot A = on M and such that the 
nonintegrable phase factor [17J for a closed path from the homotopy class g a or 
equals e 2ma or e 2m ^, respectively (assuming that < a, (3 < 1). The Hamiltonian 
then acts as the differential operator (— iV — A) 2 in L 2 (M). A unitarily equivalent 
and for our purposes more convenient possibility is to work with the Hamiltonian 
H\ = —A in the Hilbert space J^a of A-equivariant functions on M, as introduced 
in Section 12.11 Parallelly one considers the free Hamiltonian H = —A in L 2 (M). 
H is T-periodic. In order to compute, according to prescription the propagator 
}C A (t,x,y) associated to H\ one needs to derive a formula for the free propagator 
1C(t,x,y) on M. Such a formula is recalled below following [9]. 

Let "d be the Heaviside step function. For x, y G M U A U B set x{ x -> 2/) = 1 if the 
points x, y can be connected by a geodesic segment, and x( x , y) — otherwise. Given 
in addition t G R we define 

Z(t,x,y) = fi{t)x(x,y) ex P^ dist 2 (x,?/)^ , 

Furthermore, for X\, x 2 , x 3 G M U A U B such that x( x ii x 2) = x( x 2, x 3) = 1, and for 
ti, t 2 > we set 

v { x x:c )=*((*-* +^(£)T - (" - + "^fis))" 1 ) 

where 8 = Z xi, x 2 , X3 G R is the oriented angle and t\ = dist(xi, x 2 ), r 2 = dist(x2, X3). 
Note that if the inner vertex x 2 belongs to the set of extreme points AU B then the 
angle 9 can take any real value. 
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We claim that the free propagator on M equals 

fC(t, x, x ) = 2_. /C 7 (t, x, x ) (5) 

7 

where the sum runs over all piecewise geodesic curves 7 : xo — > C\ — > . . . — > C n — > x 
with the inner vertices Cj, 1 < j < n, belonging to the set of extreme points A U B. 
This means that it should hold x( x o,Ci) = x{C\,C2) = ... = x(C n ,x) = 1. Let 
us denote by I7I = n the length of the sequence (Ci,C 2 , ■ ■ ■ ,C n ). In particular, if 
I7I = then 7 designates the geodesic segment x — > x. To simplify notation we 
set everywhere where convenient Co = 20 and C n+ i = x. With this convention, the 
summands in ([5]) equal 

/C 7 (t, x, Xo) 

= [ dt n ... dt 5(t n + ... + t -t)f[v( Cj+2 .> Cj + U Cj ) fl Z(tj, C j+1 , Cj). 

(6) 

In particular, if I7I = then /C 7 (t, x, x ) = Z(t, x, x ), and if I7I = 1 then 7 designates 
a path composed of two geodesic segments x — > C — > x, with C G A U £>, and 

/C 7 (t,x,x ) =$( t )f o V ( y X t '^ , s X ° s } Z(t-s,x,C)Z(s,C,x )ds. 

3.2 Auxiliary relations 

As is well known, in R 2 it holds true that 

dx dy) x + iy d(x)5(y) 
and, consequently, 

x + iy 

With the aid of the last equality one can verify that the relation 

+ ; | + ^) (« + < "* (;))" = 7? (*(« - o*w - - 

(7) 

is valid on the domain i > 0, r > 0, 9 G R. It is straightforward to see that, on the 
same domain, 
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Combining (J7|), §E§ and the Leibniz rule one finds that 

d d 2 Id 1 d 2 \ / , ftW 1 1 / .r 2 



dt Or 2 ~^ r dr ~^ r 2 dO 2 ) \ ^~ °^ \r ) ) t 6 \' At 
= ^ exp fi Q (8(t - r)8'(6) - ir5'(t - r)5(9)) . (9) 

Equipped with ([9]) one can prove the equality 

d d 2 Id 1 d 2 \ f l / , f(t-s)r 



dt dr 2 r dr r 2 d6 2 J J \ \ sr 



2ixr$ f (r + ?"o)r 



\ ex P 1 
r / (r + ro) \ At 



t-s 



1 r + r 



r ff 1+ . r (r + r )\ J Jr±_\ + Jr^ f Jjr^ ]] m 



r + r o\\ At J \r + r / r + r \r + r 

"(10) 

which is true in the sense of distributions for any r > and / G C 1 ([0, +oo[), again 
on the domain t > 0, r > 0, 9 eR. Note that 



In particular, letting 



j exp (i ^ as e -> + in ]0, +oo[ ). 



/( a ) = ±«p(i 



one derives the following equality which is true in the sense of distributions, 
d d 2 Id 1 d 2 \ [ t f n /(^-s)r xW 



<9t <9r 2 ^ r 9r ^ r 2 <96* 2 ) Jo \ ^ ^ V sr 

x (T i i); exp (' i (3(r : ^) + s)) d5 

Let us also recall the following basic fact concerning the generalized Laplacian. If 
G C M is an open set with a piecewise smooth boundary, xg is the characteristic 
function of G, if is the normalized outer normal vector field on dG and r\ is a smooth 
function on M then, in the sense of distributions, 

dr) d 

HvXg) = (^v)Xg ~ ^S dG - —(r]5 dG ). (12) 

dn an 

The distribution 5qg is the layer supported on the curve dG which is defined by the 
curve integral 

^eC^(M), 5 dG {^)= [ ipdi. 

JdG 
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3.3 Verification of the formula 

We have to show that, for xq G M fixed, the propagator K.(t,x,Xo) defined in ([5]), (J6j) 
verifies the condition 

d \ 

i— + A\ K(t, x, x Q ) = i 5{t)5(x, x ) onKxM. 
This is equivalent to showing that 



and 



lim fC(t,x,xo) = 8(x,xo) (13) 
t — >o+ 

(i-^ + A J JC(t, x, x Q ) = for t > 0, x G M. (14) 

Equality (TTH|) is rather obvious. Since Z(t,x,x ) looks on the sheet {x; x( x -> x o) — 
1} as the free propagator on IR 2 we have 

lim Z(t,x,x ) = 5(x,x ). 

t — >0-)- 

From a similar reason, lim^ 0+ Z(t, x,C) = if C G AU B and x runs over M. Hence 

lim fC 7 (t,x,x ) = if |7| > 1. 

t >0-L 



Let us proceed to the verification of (1141) . First we introduce some notation related 
to the geometry of the universal covering space M. Denote by g the distance dist(a, b). 
Observe that if d, C 2 G A U B then X (Ci, C 2 ) = 1 if and only if dist(Ci, C 2 ) = g. If 
this is the case then necessarily G\ G A and C2 G B or vice versa. 

For ac G M U A U B set 

D(a;) = {i/GM;x(a: J 3/) = l}. (15) 

If x G M then D(x) can be identified with the aid of the projection tt : M — > M with 
the plane cut along two halflines with the limit points a and b, respectively. Thus the 
border of D(x) consists of two pairs of halflines. One pair has a common limit point 
A G A and is denoted dD(x; A), the other pair has a common limit point B G B and 
is denoted dD(x; B). We have 

dD(x) = dD(x; A) U dD(x; B). (16) 

If C G A U B then D(C) resembles the universal covering space in the one- vortex 
case. It can be viewed as a union of countably many sheets glued together in a 
staircase- like way. Each sheet contributes to the border of D(C) by a pair of halflines 
with a common limit point C . Thus the border dD(C) is formed by a countable union 
of pairs of halflines: if C G A then we write 

dD{C)= |J dD(C;C), (17) 
C"eB,dist(c,c)=e 
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if C G B then 

8D(C)= |J dD{C-C). 

C'eA,diBt(C,C')= e 

Let us first examine the case |7| = 0. It holds 

d \ 
i— + A)Z(t,x,x o ) = 0. 



dt 

for t > and x G -D(x ). Observe also that 

d 



dlt 



Z(t ,x,Xq) = for x G dD(xo) 



where it is the normalized outer normal vector field on the border dD(xo). This is so 
since, in the polar coordinates centered at Xq, Z(t, x, xq) does not depend on the angle 
variable. Let us also note that the function Z(t,x,xo) can be continued smoothly in 
the variable x over the borderline of the domain D(xq). Thus, in virtue of ffT^j) . we 
obtain (for t > 0, x G M) 

( d \ d 

[ i— + A J Z(t, x, x ) = ~Q=^ { z (t, x , x o) S 9D (x )) ■ (19) 

Remark. In (1191) as well as everywhere in this section we use the following convention. 
The value of a density (which is in this case Z(t,x,xo)) on the border dD(xo) is 
understood as the limit achieved from the interior of the domain D(xq). 

Next we discuss the case I7I = 1. Then 7 designates a piecewise geodesic curve 
x — > C — > x, with C G A U B. Denote by 7' the geodesic segment Xq — X (provided 
x G D(xq)). We have 



/C 7 (t,x,x ) = —-x(x,C)x(C,x 




(t - s)r \\ ( (it- s)r 

11 1 6 + 7T + z log 1 



sr J J \ \ sr 



(t - s)s 



exp ( z (i(TT7) + S)) ds 



(20) 



where r = dist(x,C), r = dist(C, x ) and 6 1 = Zx Q ,C, x. 

An application of the differential operator (id t + A) to the RHS of (|2"U1) in the 
sense of distributions results in several singular terms supported on one- dimensional 
submanifolds. First, due to the discontinuity of the characteristic function x(x, C), 
the application of the Laplace operator leads to two terms supported on the boundary 
dD{C) (see (fT2|) ). Second, as it follows from (fTTT) . the singularity of the integrand for 
the values 9 = ±7r and r /s = r/(t — s) produces terms supported on the submanifold 
determined by 9 = ±7r, and this set is nothing but a part of the boundary of the 
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domain D(x ), namely dD(x ; C). Notice that for 9 = ±n it holds r + r = dist(x, x ) 
and d/dlt = ±r~ 1 d/d9. Moreover, in the polar coordinates centered at C, 

$dD(x ;C) = - (5(9 - tt) + 6(9 + tt)) . 

Thus the latter contribution takes the form 

1 d ( ( i ,. „ , \ ,,.,„ x „,„ ,A d 



eXp (i dist(x, 2; ) 2 J (8(6 -tt) -5(0 + 7r)) j = (tC t ' (x, x )6 dD{x0jC )j 



Anitr 2 89 

where /Cy(t,x,x ) = Z(t,x,x ). In summary, we obtain 
(i^ + &j /C 7 (t,ar,j? ) = - /C 7 (t, x, x )^ <W(G) - ^ (£ 7 (t, x, x )<W(C)) 

+ t—; (fcy(t,x,xo)SdD{x 0) c)) ■ (21) 

Finally let us consider the case |7| > 2. Thus 7 is a piecewise geodesic curve 
x — ► Ci — >• . . . — > C n — > x, n > 2. Denote by 7' the truncated geodesic curve 
x — > Ci — ► . . . — ► C n _i — > x (provided x G D(C n _i)). One can express 

K. J (t,x,x ) = / dt n _i . . .dt V[ X l ™ C . n ~ X ) Z(t - t,x, C n )F 7 (t , . . . ,t„_i,x ) 

= ^~x(x,C n ) / dt n _ 2 ...dt / dt n _i 
^ Jr™- 1 Jo 




1 



+ 71 + Z I06 



(t-T'-t n -i)Q 

t n -ir 

t-T>- t n _ X ^ 4(t - T> - f~0 J ^ • • • ' ^' (22) 

where 

r = t + • • • + * n -2 + <n-i, r' = t + • • • + *»-2, >" = dist(C n , x), 9 = 1. C n _i, C n , x, 



and 



n-2 , s n-1 

F 7 (t ,...,t n _ 1 ,x )=n^ j i : j / j n^^+^i)- 



j= o v 7 i=o 

An application of the differential operator (i<9 t + A) to the RHS of ([2"2"]l in the sense 
of distributions again produces several singular terms. In consequence of the disconti- 
nuity of the characteristic function x(x, C n ) a single and a double layer supported on 
the boundary dD(C n ) occur (see (fl2~]l ). The singularity of the integrand for the values 
9 = ±n and 1 = r/(t — r f — t„_i) produces terms supported on the part of the 
boundary of the domain D(C n _x), namely on dD(C n -i, C n ). This time one can apply 
identity ( ITUl) . In order to treat the resulting terms the following equalities are useful. 
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Suppose that 9 = ±n and so x G dD(C n -\\ C n ). Set 

r' = r + q = dist(C n _i, x), 9' = Z C„_ 2 , Cn-i, x. 

If e/i„_i = r/(t - r' - t n _i) then 

e(t-r') t-r'-t^ t-r' 
t n -i = ; and = — . 



Moreover, 



and 



7, ex P 



ir 



V 



4(t - r) 

gs 2 /r', si 



Z (t n -i, C n , C n _i) — Z(t — r , x, C n _i) 



V 



S2, Si 



Observe also that 
9s 



exp 



4(t - r' - s) 



exp 



4s 



=g(t-T>)/r' 



exp 

and for 9 = tt, 



ir 



4(t - r' - s) 



exp 



z^r 
47 



exp 



s=g(t—T')/r' 



% r 



4(t - t') 



9 yf C n , C n _i, C n _2 

s, t n -2 



<9s 



s= e (t-T>)/ 



T , )/r , g(t-T')d9' V t-r',t n „ 2 



A similar relation holds true also for = — 7r. 

After a bit tedious but quite straightforward manipulations one arrives at the final 
equality 



d 

i— + A ) /C 7 (t,x,x ) 



d 
dlt 

d 
dlt 







1Cy(t,x,xo) ) 6 9D ( Cn ) ~ (]C r (t,x,xo)6 dD (c n )) 



JCj'(t,x,x ) ) S, 



>dD{C n -r,C n ) 



d 



+ ^ (fc 7 '(t,x,x )5 gD{ c n - i; c n )) ■ 



(23) 



Now we can show equality (1141) when taking into account ( fl9i) . ( 1211) and ( 1231) . It 
is true that 



13 



i-^ + A) !C(t..r..r {) ) 



E 

M>2 



d \ d 

JC 7 (t,x,x ) S aD{Cn ) - {fc. 7 (t,x,xo)8dD(Cn)) 



d \ d 

+ ( ^/Cy(t,a;,Xo) J <W(C n _ i; C„) + g=> {^'{t, X, X )5 dD{Cn . 1 ;C n )) 



d \ d 

JC^(t,X,X ) 5 d D{G) - (K--y(t, x , x o)8dD(c)) 



n J v n 



(9 n : 
= 

where we have used (fTBj) . (JT71) and (fTH|) . 



(^(t,ar,ar )<yar>(*o)) 



4 The Aharonov-Bohm effect with two vortices: 
the propagator associated to H\ 

Without loss of generality we can suppose that the vortices are located in the points 
a = (0, 0) and b = (g, 0). Let (r , 6 a ) be the polar coordinates centered at the point a 
and (r&, 0&) be the polar coordinates centered at the point b. To express the propagator 
for it is convenient to pass to a unitarily equivalent formulation. Let us cut the 
plane along two half-lines, 

L a = ] — oo, 0[ x{0} and L b = }g, +oo[ x{0}. 

The values 9 a = ±7r correspond to the two sides of the cut L a , and similarly for dj, and 
Li. The unitarily equivalent Hamiltonian H' A is formally equal to —A in L 2 (IR 2 ,d 2 :c) 
and is determined by the boundary conditions along the cut, 

^(r«, 9 a = tt) = e 2 ^>(r a , 6 a = -tt), d r >(r a , # a = tt) = e 2 ™<9 ra V(ra, #a = -tt) , 
^(r 6 , ^ = tt) = e 2 ^ ^ = -tt), <9 r ,^, 9 h = tt) = e 2 * ip d rb ij)(r h , fl 6 = -tt) . 

In addition, one should impose a boundary condition at the vortices, namely ip(a) = 
if)(b) = 0. 

Let us denote D = R 2 \ (L a U L b ). Then one can embed D G M fundamental 
domain. We wish to find a formula for the propagator /C' A (t, x, x ) associated to the 
Hamiltonian H'\. It can be simply obtained as the restriction to D of the propaga- 
tor K, A (t,x, xq) associated to the Hamiltonian H\. On the other hand, to construct 
JC A (t,x, xq) one can apply formula ([3]) and the knowledge of the free propagator on 
M, see (ED, (E]). Thus we get 

/C A (t, x, x ) = J2J2 M*» 5 • a;, *o). (24) 

9 er 7 
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Fix t > and x , x G -D. One can classify piecewise geodesic paths in M, 

7 : x -»• Cx -»• . . . -»• C„ -»• </ • x, (25) 

with Cj G .A U £> and g G T, according to their projections to M. Let 7 be a finite 
alternating sequence of points a and b, i.e., 7 = (ci, . . . , c n ), Cj G {a, 6} and Cj 7^ Cj + \. 
The empty sequence 7 = () is admissible. Relate to 7 a piecewise geodesic path in 
M, namely Xq — ► c\ — > . . . —> c n —>■ x. Suppose that this path is covered by a path 
7 in M, as given in (|25p . Then Cj G A iff Cj = a and Cj G £> iff Cj = b. Denote 
the angles Zx ,ci,c 2 = and Zc n _i,c n ,x = 6. Then the angles in the path 7 
in ( |25i) take the values Zxq,Ci,C2 = + 2irki, ZC n -i,C n , g ■ x = 9 + 2iik n and 
Z Cj, Cj + 2 = 2nkj + i for f < j < n — 2 (if n > 3), where ki,...,k n are integers. 
Any values k\, . . . , k n G Z are possible. In that case the representation A applied to 
the group element g occurring in (1231) takes the value 

A(g) = exp(2m(k 1 a 1 + . . . + k n a n )) 

where Oj G {a, /3} and Oj = a if Cj = a, and cr,- = (3 if Cj = 6. 
Using the equalities 

^2 ex v( 2mak ) [q + - 7T + i s ~ ^ + 2vrA; + vr + is J 
fcgZ v 1 1 / 

/- +oc exp((fl + i g )r) ^ 
= -siii to) / — — — — — -dr 

J_ 00 sm(vr(a + irj) 

and 

r exp((fl + igjr) d _ 2 exp(-a(s - ifl)) 
y.oo sin(7r(a + ir)) 1 + exp(-s + \9) ' 

that are valid for < a < 1, |0| < 71", one can carry out a partial summation in (J24l) 
over the integers ki, . . . , k n . This way the double sum in (EMI) reduces to a sum over 
finite alternating sequences 7. 

Let us conclude our contribution by giving the resulting formula for /C /A (t, x, x ). 
We set 

C a = 1 or (a = e 27ria or ( a = e" 27riQ 

depending on whether the segment x$x does not intersect L a , or xqx intersects L a and 
Xo lies in the lower half-plane, or xqx intersects L a and xq lies in the upper half-plane. 
Analogously, 

£ 6 = 1 or C b = e 2 " i/3 or & = e~ 2nil3 

depending on whether the segment ~XqT does not intersect L b , or xqX intersects L b and 
xo lies in the upper half-plane, or x^x intersects L& and x lies in the lower half-plane. 
Furthermore, let us set 

(a = e iar >% Cb = e i ^\ where 77 , 776 G {0, 2 vr, -2 vr}. 
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Then one has 



K,' A (t,x,x ) 

j. /■ 1 ( .\X-XQ 

= CaCb -7—1 eX P 1 



Anit \ At 

CO 

An 2 i Jo ti Jo t 
r 

Ati Ato) ) 1 + exp(-s a + i9 a - i9 0a ) 
U 

Jo ti Jo to 



y I ; - jj_ + r L\\ exp[-a(s a - i(9 a - 6 0a - Va)] 



x exp \i\ -r- + 



y 2 , vi 2 

'6 , ' 06 



exp[-p(s b - i(9 b - 6ob - Vb)] 



Ati Ato) J 1 + exp(-s b + i9 b - i9 0b ) 



7,n>2 



x exp(l(^ + -^ + ... + f + f) ) S^.0.f) (] ). 



where 



sin(7ra w ) exp[-a n (s ra - id)] sin7rcr n _i exp(-a ra „is n _i) 
7r 1 + exp(-s n + id) it l + exp(-s„) 
sin(7ra"i) exp[— <7i(si — i9 )] 



x . . . x 



7r 1 + exp(-si + i6 ) ' 



and 



s a = log( t -f^ L j , s b = log( ) , Sj = log( j for 1 < j < n. 

\t r a J \t r b J Vj-iTjJ 

In addition, (r, 9) are the polar coordinates of the point x with respect to the center 
c n , (r ,9o) are the polar coordinates of the point x with respect to the center c\. 
The sum n>2 runs over all finite alternating sequences of length at least two, 
7 = (ci, ...,c n ), such that for all j, Cj G {a, b}, Cj ^ c?+i, and <jj = a (resp. f3) 
depending on whether Cj = a (resp. b). 
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